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ABSTRACT 

Recently, Kallivayalil et al. derived new values of the proper motion for the Large and Small 
Magellanic Clouds (LMC and SMC, respectively). The spatial velocities of both Clouds are 
unexpectedly higher than their previous values resulting from agreement between the available 
theoretical models of the Magellanic System and the observations of neutral hydrogen (HI) 
associated with the LMC and the SMC. Such proper motion estimates are likely to be at odds 
with the scenarios for creation of the large-scale structures in the Magellanic System suggested so 
far. We investigated this hypothesis for the pure tidal models, as they were the first ones devised 
to explain the evolution of the Magellanic System, and the tidal stripping is intrinsically involved 
in every model assuming the gravitational interaction. The parameter space for the Milky Way 
(MW)-LMC-SMC interaction was analyzed by a robust search algorithm (genetic algorithm) 
combined with a fast restricted N-body model of the interaction. Our method extended the 
known variety of evolutionary scenarios satisfying the observed kinematics and morphology of the 
Magellanic large-scale structures. Nevertheless, assuming the tidal interaction, no satisfactory 
reproduction of the H I data available for the Magellanic Clouds was achieved with the new proper 
motions. We conclude that for the proper motion data by Kallivayalil et al., within their 1-er 
errors, the dynamical evolution of the Magellanic System with the currently accepted total mass 
of the MW cannot be explained in the framework of pure tidal models. The optimal value for 
the western component of the LMC proper motion was found to be /Lt^jL > 
of tidal models. It corresponds to the reduction of the Kallivayalil et al. value for /i^ 
in its magnitude. 
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1. Introduction 

The discussion of the origin and evolution of 
the Magellanic System has become very intense 
since the new proper motion data for the selected 
LMC/SMC stars were acquired by the Hubble 
Space Telescope (HST). The HST measurements 
yielded the new values for the mean proper mo- 
tions of the Magellan ic Clouds with an unpre ce- 
dented accuracy (see iKallivavalil et all l2006al lbh. 



In compari son with the previous observational 
studies (e.g.lJones et alj|l994t iKrouoa et al.lll994 ; 
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Kroupafc BastiarJ " 1997 ). the corresponding mea 
surement errors were reduced by a factor of 10. 
Even though the latest proper motion estimates 
are consistent with the previous observational re- 
sults within the 1-a errors, their actual position 
in the velocity space of the Clouds is quite unex- 
pected. 

The current velocities of the LMC and the 
SMC are critical input parameters of any evo- 
lutionary model of the System. However, regard- 
ing the large heliocentric distance to the Magel- 
lanic Clouds, nobody attempted the observational 
determination of their proper motion until the 
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paper s by I Jones et"aL ( 1994 ) and Kroupa et al 
(1994) were carried out. Some more studies have 
contr i buted to the research (e.g.lKroupafeBastianl 
19971: brake et al.ll200ll IPedreros et alj l2002h of- 
fering a span of the mean proper motion val- 
ues, but reaching no substantial improvement 
in the measurement precision, that was still of 
the same order as the derived values themselves. 
Such errors adm it a wide variety of scenarios for 
the interaction ( Ruzicka et al. 2007h . and make 
the observational estimates serve only as quite 
weak tests of the results based on the theoretical 
studies of the MW-LMC-SMC interaction (e.g. 
LinfcLvnden-Belllll98ilGardiner et al.lll994h . 

The theorists focusing on the evolutionary his- 
tory of the Magellanic System have widely agreed 
on two basic physical processes dominating the 
formation of the Magellanic large-scale struc- 
tures, including the Magellanic Strea m, the Bridge 



and t he Leading Arm (for details see lBriins et al 
2005). Those are the tidal fields and the ram 
pressure stripping. 

The tidal origin of the ext ended Magellanic 
struct ures was investigated by iFujimoto k Sofue 
(|l976t) . who assumed the LMC and the SMC to 
form a pair gravitationally bound for several Gyr, 
moving in a flattened MW halo. They identi- 
fied some LMC and SMC orbit al paths leading to 



the creation of a tidal tail. Lin k Lvnden-BeU 



(|1977l ) pointed out the problem of the large 
parameter space of the MW-LMC-SMC in- 
teraction. To reduce the size of the param- 
eter space, they neglected both the SMC in- 
fluence on the System and dynamical friction 
within the MW halo, and showed that such 
a configuration allows for the existence of a 
LMC trailing tidal s t ream. Following studies 
bv iMuraifcFuiimotol (Il980h. iMuraifcFuiimoto 
' l984h - iLin k Lvnden-Belll (Il982h . iGardiner etal 



19941 ). or ILin et al.1 ( 19951 ) extended and devel- 



oped tidal models of the MW-LMC-SMC inter- 
action. Generally speaking, the tidal mechanism 
becomes efficient enough if the timescale for the 
interaction is several Gy r ( Gardiner et al. 19941 
Gardiner k Noguchilll996l ). 

Meurer et al. (jl985h involved continuous ram 



pressure stripping into their simulation of the 
Magella nic Sy s tem This approach was followe d 
later bv lSofuel (| 19941 ) or bv lMoore k David (|l994l ). 
who simplified the interaction between the LMC 



and the SMC, however. The Magellanic Stream 
was formed of the gas stripped from outer re- 
gions of the Clouds due to collisions wi t h the 
MW extended ionized disk. iHeller k Rohlfel (| 19941 ) 
argue for a LMC-SMC collision resulting into 
the gas distribution to the inter-cloud region 
where it was stripped off by ram pressure as the 
Clouds moved throug h the hot MW halo. Re- 
cently, iBekkifcChibal (|2005l ) applied a complex 
gas-dynamical model including star-formation to 
investigate t he dynamical and chemic al evolution 
of the LMC. lMastropietro et all (|2005l ) introduced 
their model of the Magellanic System including hy- 
drodynamics (SPH) and a full N-body description 
of gravity. They studied the inter action between 
the L MC and the MW. Neit h er iBekki fc Chiba 
( 2005 ) nor iMastropietro et ah ( 20051 ) considered 
the SMC gas in their models. However, it was 
shown that the Stream, which sufficiently repro- 
duces the observed H I column density distribu- 
tion, might have been created without the SMC 
gaseous compone nt or even without the LMC- 
SMC interaction (jMastropietro et al.l[2005h . The 
history of the Leading Arm was not investigated. 

In general, hydrodynamical models allow for a 
better reproduction of the H I column density pro- 
file of the Magellanic Stream than tidal schemes. 
However, they constantly fail to reproduce the 
Magellanic Stream radial velocity measurements 
and especially the high negative velocity tip of 
the Magellanic Stream. Both families of mod- 
els suffer from serious difficulties when modeling 
the Leading Arm. Similar requirements as for the 
tidal models hold for the ram pressure stripping 
schemes concerning interaction timescales, unless 
the density of the extended gaseous halo of the 
MW is increased substantially over its observa- 
tional estimates, amplifying the hydrodynamical 
interaction. The ram pressure force is propor- 
tional to v 2 (relative velocity of the interacting 
gaseous objects). Thus, the process of gas strip- 
ping becomes more efficient as the velocity in- 
creases. However, for finite-size objects, it is ac- 
companied by shorter interaction timescales due 
to the reduced crossing time. 

Here we came to the actual point of contro- 



versy r elated to the papers by iKallivavalil et al 



(|2006al lbT): the HST proper motion values put the 
Clouds on highly eccentri c (or even unbound ) or- 
bits around the Galaxy. Besla et "aD ( 20071 ) an- 
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alyzed the orbital motion of the LMC using the 
HST data and brought convincing arguments jus- 
tifying the previous statement. Such a result may 
have serious consequences for the proposed forma- 
tion mechanisms of the Magellanic System, since it 
strongly discriminates the tidal scenario and prob- 
ably also the ram pressure-based models. To ex- 
plain the evolution of the Magellanic System by 
either of the mentioned processes, the negative 
total energy of the Clouds on their orbits about 
the MW is needed, which corresponds to multi- 
ple perigalactic approaches over the Hubble time. 
It is highly desirable to verify the reliability of 
the available models of the MW-LMC-SMC in- 
teraction if the orbital angular momentum of the 
Clouds is as high as found by iKallivavalil et al" 
(|2006alfbh . 



Our paper presents the results of the search 
of the parameter space for the MW-LMC-SMC 
interaction dominated by tides. The approach 
applied was introduced by iRuzicka et ail (|2007l ) 
and is based on an evolutionary optimization of 
the model input according to its ability to re- 
produce the HI observati ons of the Magellanic 
System (IBruns et al.l 120051 ). The predictions by 



Besla et all (j2007l) and by many others regard- 
ing the ins ufficient performance of t idal models in 



case of the IKallivavalil et alJ (|2006aU bl) proper mo- 



tions can only be confirmed if the entire parameter 
space for the interaction is explored. The method 
of the automated search for good models by the ge- 
netic algorithms (GA) enabled us to perform such 
an analysis for the first time. In addition to the 
LMC/SMC velocity problem, we also intend to 
answer two more questions: Are all the studied 
parameters of the same importance for successful 
modeling the MW-LMC-SMC interaction? Does 
the evolution of the System show similar behavior 
over various scales in the parameter space? 

2. Parameter space of the interaction 

Our model is built in a galactocentric Carte- 
sian frame, assuming the present position of the 
Sun r Q = (— 8.5, 0, 0) kpc and i ts spatial velocity 
v m = (10.0,225.2, 7.2) kms" 1 (jPehnen fc Binnev 
1998). In total, our study involves over 20 inde- 
pendent parameters, including the initial condi- 
tions of the LMC and the SMC motion, their to- 
tal masses, parameters of mass distribution, parti- 



cle disk radii, and orientation angles, and also the 
MW dark matter halo flattening parameter. Some 
of the parameters were constrained by theoretical 
studies (including scale radii e of the LMC/SMC 
halos, the Coulomb logarithm A for the dynamical 
friction in the MW halo, and the halo flattening 
parameter q for the model of the MW gravitational 
potential). Their mean values and sear ched errors 
were discussed in IRuzicka et al.l ( 2007 ). This sec- 
tion focuses on the observationally estimated pa- 
rameters of the MW-LMC-SMC interaction. 

Figure [1] shows the current proper motion of 
the LMC and the SMC as estimated by var- 
ious observational methods, together with the 
portion of the velocity space that we studied. 
For convenience, the proper motion vectors were 
decomposed into the norther n (m n ) and th e 



wester n (ji ) components (see IKallivavalil et al 



2006af ). We included every case that is accept- 
able up to date. In case of t he LMC pro p er mo 



tion, the 



Kroup a et al 



mea surem ents 
(Il994h 



by I Jones et al. 



Kroupa & Bastia 



Drake et alJ 



l200lh: iPedreros et al 
2006 



i 



(2002) 



1994); 



19971) : 



and 
The 



IKallivavalil et al.l (|2006al ) were considered 
proper m otion errors for the SMC a r e based on th e 
results bvlKroupafcBastianl (Il997l):llrwinl (Il99 



d2003l): [Anderson fc Kind (|2004all 
and IKallivavalil et al.l (|2006rJ ) 



iFreire et al 



Efficiency and reliability of the optimization 
method - genetic algorithms - always depend on 
the number of possible solutions, that is extremely 
high in case of the Magellanic parameter space. 
To resolve the difficulty we performed a two-level 
search. First, the full volume of the parameter 
space was explored (referred to as a "global scale" 
hereafter) to obtain a low resolution information 
about the behavior of the tidal model for vari- 
ous combinations of the input parameters. Sub- 
sequently, we reduced the volume of the param- 
eter space by a factor of 10 4 ("local scale" here- 
aft er) to study only the 1-a proper motion errors 
by IKallivavalil et all (|2006alibh with high resolu- 
tion. Finally, both searches were compared. The 
extended ranges of the LMC and the SMC proper 
motion were 
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(-0.50, +1.80) mas yr" 1 (1) 

(-2.23, -0.70) mas yr" 1 (2) 

(-2.50, -0.46) mas yr^ 1 (3) 

(-2.05, +0.00) mas yr" 1 . (4) 
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To perform the high-resolution search, the vol- 
ume of the parameter space was reduced by 
adopting the proper m otion values derived by 



adopting tne proper rn otic 
iKallivavalil et alj (|2006al lbh: 



rime 
Mime 

rsmc 

w 

rsmc 



(+0.39, +0.49) mas yr" 1 
(-2.11, -1.95) mas yr^ 1 
(-1.35,-0.99) mas yr" 1 
(-1.34,-0.98) mas yr" 1 . 



(5) 
(G) 
(7) 
(8) 



The above introduced proper motion space for the 
Magellanic Clouds is illustrated by Fig. [TJ 

Unlike the proper motion of the Clouds, their 
LSR radial velocities could be measured w i th hig h 
accuracy. Following Ivan der Marel et al.l ([2002) , 
we set v^ c = 262.2±3.4kms- 1 The SMC radial 
veloci ty error was estimated by lHarris fe Zaritskyl 



(l2006h as i>££ = 146.0 ± 0.6 km s" 1 . 

The helioc entric position vector of the L MC was 



adopted from Ivan der Marel et al 



( 20021) . i.e. the 
S lmc ) = (81.90° ± 
modulus is 



equatorial coordinates are (aq mc , 
0.98°, -69.87° ± 0.41°), its distance 
(to - M) lmc = 18.5 ± 0.1. The equatorial co- 
ordinates of the SMC were set to the ranges 
'- w,<Uc) = (13 .2° ± 0.3°, -72.5° ± 0.3°) (see 



Stanimirovic et al.ll2004 and references therein). 
Van den Bergh ~ (f2000l ) provided a great compila- 
tion of various distance determinations for the 
SMC, and we used his resulting distance modu- 
lus (to - M) smc = 18.85 ± 0.10. 

Several observational determinations of the 
LMC disk plane orientation hav e been published 
so far (see, e.g. Lin et al. 19951 ). In our param- 
eter study, the LMC inclination i and position 
angle p together with th eir errors agree with 



van der Marel et al 



( 20021) . i.e. % = 34.7° ± 6.2° 
129.9° ± 6.0°. As the SMC misses a 
disk, the orientation and the po- 
to the SMC "bar" 
Based 



and p 

well defined 

sition angl e usually refer 
defined by iGardiner fc Noguchil (1 19961) . _ 
on the estimates by IVan den Berghl ( 2000h or 



Stanimirovic et al 



we adopted the error 
ranges % = 60° ± 20° and p = 45° ± 20° for the 
SMC initial disk inclination and position angle, 
respectively. 



Gardiner et alj (|1994T ) analyzed the H I surface 



contour map of the Clouds to estimate the initial 



LMC and SMC disk radii 



r disk 
lmc 



and 



„disk 



respec- 



tively. Regarding the absence of a clearly defined 
disk of the SMC and possible significant mass re- 



distribution in the Clouds during their evolution, 
the results require a careful treatment. 
Current to tal masses tow and 



aim // i. S mc follow the 
estimates by IVan den Berghl (2000). The masses 
of the Clouds are functions of time and evolve 
due to the LMC-SMC exchange of matter, and 
as a consequence of the interaction between the 
Clouds and the MW. Our test-particle model does 
not allow for a reasonable treatment of a time- 
dependent mass-loss. Therefore, the masses of the 
Clouds are considered constant in time, and their 
initial values at the starting epoch of simulations 
are approximated by the current LMC and SMC 
masses. 

The dynamics of the MW-LMC-SMC interac- 
tion is critically dependent on the density distri- 
bution and the total mass of the Galaxy. We 
model the MW by the simple axially symmet- 
ric logarithmic potential involving 3 parameters 
(jBinnev fc Tremainel [l987h . Only the MW halo 
flattening q was treated as a free parameter in 
this stu dy, varying within t he range (0.78, 1.20) 
(see also lRuzicka et al.l 120071 ). and thus introduc- 
ing a spread in the total mass of the Galaxy 
m M w{q) = (0.92, 2.15) • 1O 12 M within the radius 
of 200 kpc. 

3. Methodology 

We investigate the pure tidal models in this 
paper, as they were the first ones devised to 
explain the creation o f the Magellanic Stream 
(|FuiimotofcSofuelll976h . Reside that, the tidal 
stripping is intrinsically involved in every model 
assuming the gravitational interaction. The model 
itself is an adva n ced y ersion of the scheme by 
FuiimotofcSofuel (|l976l) : it is a restricted N- 
body (i.e. test particle) code describing the grav- 
itational interaction between the Galaxy and its 
dwarf companions. The potential of the MW is 
dominated by the flattened dark matter halo, and 
the dynamical friction is exerted on the Magellanic 
Clouds as they move through the halo. The LMC 
and the SMC are represented by Plummer spheres, 
initially surround ed by test-particle d isks. For 
further details see Ruzicka et al.1 ( 2007 ). 



3.1. Genetic algorithm search 

As mentioned already, the search itself was per- 
formed by GAs that mimic the selection strategy 
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of the natural evolution. Holland first pro- 

posed the application of such an approach on opti- 
mization problems in mathematics. Recently, per- 
formance of GAs was stud i ed for galax ies in in- 



teraction (sec 



Wahddll998t iTheid Il999h. As an 



example, TheisfeKohld (|2001 ) analyzed the pa- 



rameter space of two observed interacting galax- 
ies - NGC4449 and DDO 125. GAs turned out 
to be very robust tools for such a task if the 
routine comparing the observational and modeled 
data is appropriately defined. The approach by 



Theis fc Kohld (|200ll ) was later adopted and im- 



proved in order to explore the inte raction of the 



Mage llanic Clouds and the Galaxy (Ru zicka et al 
l2007h . The comparison between the model and 
observations became more efficient by involving 
an explicit search for the structural shapes in the 
data. Also the significant system-specific features 
(such as a special geometry and kinematics) were 
taken into account, further improving the GA per- 
formance for exploration of the MW-LMC-SMC 
interaction. More detailed information is to be 
found in Sec. 13.21 of this paper. 

3.2. Comparison between model and ob- 
servations: the Fitness Function 



The proposed automatic search of the parame- 
ter space is driven by the routine comparing the 
modeled and ob served H I distribu tion in the Mag- 
ellanic System ( Brims et al.ll2005h . The match is 
measured by the fitness function (F) which is, in 
fact, a function of all input parameters, since every 
parameter set determines the resulting simulated 
HI data-cube. The devised function F returns a 
floating-point number between 0.0 (complete dis- 
agreement) and 1.0 (perfect match), and consists 
of three different comparisons, including search for 
structures and analysis of local kinematics. 

Efficiency of the GA is c ritically dependent on 
the applied fitness function. lTheisfeKohiel(|200ll) 
proposed a generally applicable technique based 
on comparing the relative intensities of the corre- 
sponding pixels in the modeled and observed data- 
cubes. Such a fitness function was successfully 
use d to analyze an intera ction involving two galax- 
ies ([Theis fc Kohldl200l[ ) . The mentioned compar- 
ison scheme became one of the three components 
of the fitness function developed for this project. 
Both modeled and observed H I column density 
values are scaled relative to their maxima to in- 



troduce dimensionless quantities. Then, we get 



Fi 
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where cr° b fc s , 



ijk 



(9) 

are normalized column densi- 
ties measured at the position [j, k] of the i-th ve- 
locity channel of the observed and modeled data, 
respectively. N v = 32 is the number of sep- 
arate LSR radial velocity channels in our data. 
(N x ■ N y ) = (32 • 64) is the total number of posi- 
tions on the plane of sky for which observed and 
modeled H I column density values are available. 



Ruzicka et al.l (|2007l ) further tested the perfor- 



mance of GA for the problem of galactic interac- 
tions, and an additional comparison dealing with 
the whole data-cubes was devised. It combines 
the enhancement of structures in the data by their 
Fourier filtering with the subsequent check for 
cmpty/non empty pixels in both data-cubes. The 
corresponding component of the fitness function is 
defined as follows: 



Fo 



e e e !><;;- • \>Kr 

i=ij=i k=i 



N v Ny Ny 



max 



E E E pix-^ s , E E E pix™f 



=ii=ifc=i 



-.ij=ik=i 



(10) 

where pix° b fc s G {0, 1} and pix™£ d G {0, 1} indicate 



whether there is matter detected at the position [i, 
j, k] of the 3D data on the observed and modeled 
Magellanic System, respectively. 

Effectively, such a comparison is a measure for 
the agreement of the structural shape in the data. 
No attention is paid to specific HI column den- 
sity values here. We only test whether both mod- 
eled and observed emission is prese nt at the same 



pixel of the position velocity space. Ru zicka et al 



(l2007h showed that the search for structures signif- 
icantly improves the GA performance if the struc- 
tures of interest occupy only a small fraction of 
the system's entire data-cube (< 10% in the case 
of the Magellanic Stream and the Leading Arm) . 

iRuzicka et al. ( 2007 ) also recommended and 
successfully applied a system-specific comparison. 
In case of the Magellanic Clouds, the very typical 
linear radial velocity profile of the Stream includ- 
ing its high negative velocity tip was considered 
important. The slope of the LSR radial veloc- 
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ity function is a very specific feature, especially 
strongly dependent on the features of the orbital 
motion of the Clouds. Then, the third F compo- 
nent is defined as 

F 3 = J^z*-, (11) 

1 + 



.mod 
min 



where and 
served LSR radial velocity profile of the Magellanic 
Stream and its model, respectively. The result- 
ing fitness function F combines the above defined 
components in the following way: 



F — F1F2F3 



(12) 



In principle, the GA is able to find the global 
maximum of F (i.e. the best model over the stud- 
ied parameter space), but such a process may be 
very time-consumi ng due to the possibly slow con - 
vergence of F (see lHollandlll97a lGoldberd[l989l ). 
In order to overcome such a difficulty, we searched 
the parameter space repeatedly in a fixed number 
of optimization steps, i.e. generations of models, 
and collected 120 high-quality models for either 
global or local scale of the parameter space. Dis- 
tribution of the 120 local peaks of F helps to map 
the fitness function landscape but it does not al- 
low for conclusions on the behavior of F either 
outside or inside the volume populated by the lo- 
calized models. Therefore, in the following para- 
graphs we intend to devise a reasonable method 
for the further analysis of the fitness function F. 

At this point, the reader might ask why there 
is such an attention paid to the properties of the 
fitness function itself if it, in fact, does not seem to 
provide any physical information about the inter- 
acting system of the Galaxy and the Magellanic 
Clouds. Indeed, the function F serves primar- 
ily as a driver to the GA engine. However, the 
search for good models of the observed Magellanic 
System is efficient only if relevant astrophysical 
data are supplied as the input to F. As already 
mentioned, our study deals with detailed morpho- 
logical an d kinematic inform ation from the 21cm 
survey by iBriins et al. ( 2005 ) and with the corre- 
sponding modeled data. The fitness function then 
makes a link between the observable data and the 
initial state of the Magellanic System. Here we 
came to the benefits of spending time on studying 
the function F: in principle, it allows for identifi- 
cation of all points/regions in the parameter space 



leading to reproduction of the observational data, 
and also an appropriate analysis of its behavior 
may evaluate sensitivity of the System to varia- 
tions in different parameters, i.e. their importance 
to the evolution of the System. 

The parameter space of the interacting Magel- 
lanic Clouds is very extended. Acquiring helpful 
information about the function F is quite a de- 
manding task, but we consider it feasible once the 
goals of such an analysis are properly defined. The 
investigation of the fitness function F is supposed 
to help us to answer two questions already raised 
in Section [TJ Are all the studied parameters of the 
same importance for the properties of Fl Does F 
of the system show similar behavior over various 
scales in the parameter space? 

3.3. Application of the Fitness Function 

As the first step towards better understanding 
of the fitness function, we studied the 1 D projec- 
tions of F to the plane of the j-th parameter 



F'^^fipi,. 



•iPn)i 



(13) 



where p\, ...,Pj_i, p}+x> ••• iPn are l ne specific 
values of the parameters corresponding to the i-th 
GA fit (point in the parameter space) and the pa- 
rameter pj is varied within a given range. To quan- 
tify the sensitivity of F l (pj ) to changes in different 
variables (parameters) several functions were de- 
fined. First, we have 

(14) 



which is the relative deviation of the 1 D projected 
fitness F l (pj) from its mean value F l (pj) on the 
studied interval, where 



\ 



JVp-1 



jL jr; (F*(Pi, fc )-F*(Pi)) a (15) 



fe=0 



is the corresponding deviation from F l (jpj) ex- 
pressed for Np points pj t k of the projected func- 
tion F l (pj). We also calculate the relative change 
in F 1 (pj) as 



A 



_ 1 vaax.(F i (p j )) - min(F' L ( Pj )) 



(16) 
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The behavior of F is analyzed here in terms of 
its deviation from the reference levels, which are 
established by the mean values of the 1 D projec- 
tions F z (pj). We found such an approach partic- 
ularly useful if one wishes to distinguish between 
the large-scale and localized significant changes in 
F. 

Low values of both, A J and A J F , indicate pres- 
ence of a global plateau of F z (pj), and thus weak 
dependence of the system on the j-th parame- 
ter. If the corresponding A* F is of a significantly 
higher value, local peaks (or wells) exist. Simi- 
larly, the overall considerable evolution of F l {pj) 
is revealed by high values of both A* and A* F . In 
such a case additional information is provided by 
the ratio of the functions (fT4|) and (fl6|) . As the 
ratio approaches unity, abrupt changes in F l (pj) 
are favored over its smooth and slow evolution. 

The intention of this section is to apply the 
method for the fitness function analysis introduced 
in Sec. 13.21 to the Magellanic System. The func- 
tions Aj and A* F are helpful if features and be- 
havior of the fitness function F are studied in the 
neighborhood of an arbitrary point in the param- 
eter space. However, to answer the above raised 
questions, an approach somewhat less detailed is 
sufficient. We suggest to treat the functions A* 
and Aj F statistically and to calculate their mean 
values Aj = A* and A^p 



A) F over all 120 



GA fits for both global and local scale cases. The 
function F may be characterized by significantly 
different values of AJ or AJ F depending on the se- 
lected point in the parameter space. But for now, 
we are particularly interested in general trends in 
the behavior of the fitness function (i.e. the behav- 
ior of our model for the interaction) that should 
be expected if one studies the impact of variations 
in a selected parameter. As we will learn later, 
the identified GA fits cover a large fraction of the 
total parameter space volume. This fact also jus- 
tifies the proposed statistical treatment of the 1 D 
projections of the function F in case the results 
apply on the entire parameter space. 

4. Results 

Regarding the immense difficulties accompany- 
ing the observational measurement of the proper 
motions in case of the Magellanic Clouds, the mod- 
els of the MW-LMC-SMC interaction were used 



to d raw conclusions on the motion of the Clouds 
fe.g.lFunmotofeSofudll976UGardiner et alj|l994 
Heller fc Rohlfd 11994 [Gardiner fc Noguchil Il996) . 



The mentioned models preferred either the tidal 
or hydrodynamical interactions as the processes 
dominating the evolution of the Magellanic Sys- 
tem. Generally speaking, the proposed forma- 
tion mechanisms are not efficient enough un- 
less the Clouds orbit around the Galaxy. Sev- 
eral perigalactic approaches of t he Clouds are 



expected by the tidal m odels ((Gardiner et al 



1994t IConnors et alj [2005). Shorter timescales 
for the interaction may be sufficient within 
the ram pressure scenar io dHellerfc Rohlfd Il994t 



Mastropietro et alj 120051) . However, the proper 
motions by Kallivavalil et al.l ( 2006a. b) lead to 
timescales further dramatically reduced, as the 
Clouds sh ould be approachi ng the Galaxy for the 
first time ( Besla et al.ll2007t ). The research of the 
dynamical evolution of the Magellanic System is 
at the point where our theoretical understanding 
of the MW-LMC-SMC interaction is at odds with 
some critical observational constraints. 



Ruzicka et al 



(|2007f ) realized that the previ- 
ous attempts to model the Magellanic System al- 
ways reduced the parameter and initial condition 
space for the interaction by additional assump- 
tions, such as omitting the SMC, or adopting a 
special orbit for the LMC. However, the unique- 
ness of the models is unclear unless a systematic 
analysis of the entire parameter space compati- 
ble with the available observatio ns is performed. 
The idea by iRuzicka et alj (|2007l ) was justified as 
they used a tidal model of the interaction and re- 
produced the observed HI structures for remark- 
ably different histories of the System. Such a com- 
plex approach requir es a powerful search method. 
IRuzicka et al. ( 2007 ) resolved the difficulty by em- 
ploying GAs as optimization tools characterized 
by reliability and low sensitivity to local extremes 
([Theis fc Kohldl200ll ). We adopted their approach 
to analyze the performance of pure tidal models in 
case of the Magellanic System a ssuming the LMC 
and the SMC proper motions bv lKallivavalil et al 



(|2006al lbl). 



4.1. Parameter Dependence 

The values of Aj and A^f were calculated for 
the fitness function F over both global and local 
scales of the parameter space. Table Q] presents 
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the result in the descending order according to the 
value of the function Aj. Following the previous 
explanation it is clear, that the order reveals the 
sensitivity of the fitness function to various param- 
eters. At this point, we feel qualified to answer the 
first question raised in Sec. 13.21 the parameters are 
not of the same importance for the properties of F. 
Table Q] indicates that the sensitivity of the Mag- 
ellanic System to the choice for the LMC/SMC 
proper motions is significantly higher than in case 
of any of the remaining parameters This conclu- 
sion is well illustrated by Figure[3]depicting the fit- 
ness landscape projected to the plane of the LMC 
proper motion components. The projection was 
made at the positions of the best GA fits found 
for the global and local scale cases, respectively. 

Our analysis has shown the critical dependence 
of the evolution of the Magellanic System on the 
LMC and the SMC spatial velocities. The func- 
tion A^p is not only useful for the detailed analysis 
of the fitness function, but can also serve to esti- 
mate the typical overall change in the fitness value 
over the entire studied parameter range. The value 
of A^p exceeds 0.15 for every LMC/SMC proper 
motion component on the global scale (Tab. [!}. 
Regarding the definition by Eq. (|16p . the men- 
tioned values multiplied by the factor of 2.0 tell 
us how much the ID fitness projections F l (pj) 
change typically compared to their mean values. 
As the usual mean value F l (pj) equals to ~ 0.3, 
one may assume the global proper motion ranges 
given by Eq. to ((4]) to derive the following es- 
timated changes in the fitness of our models per 
the unit step in the proper motion value: 



AFQCcVA/zIL 


~ 0.04 mas" 


- x yr 


(17) 


AF(/0/A/C c 


~ 0.06 mas" 


- x yr 


(18) 


AF( Ms N mc )/A Ms N mc 


~ 0.09 mas" 


-iyr 


(19) 


AF(vZc)/^Z c 


~ 0.06 mas" 


- x yr 


(20) 



The above listed values should be treated very 
carefully, as they simplify the real behavior of the 
fitness function. The fact that the proper motion 
of the SMC is even more critical than that of the 
LMC indicates, that both Clouds serve as sources 
of matter for the Magellanic large scale structures, 
but the SMC contribution responds to the choice 
for the orbit more strongly. It is another nice illus- 
tration of how much information about the phys- 
ical properties of the Magellanic System can be 



actually obtained from the fitness function. 

The estimates given by the Eq. (17]) to (20]) give 
a hint at the relation between the fitness and the 
physical properties of our models. Later in Sec. 14.31 
the lower limit for the fitness of the satisfactory 
models will be established and discussed. While 
the best model we found yields F best = 0.514, the 
fitness threshold level is Fu m > 0.434, i.e. all 
the acceptable models are found within the fit- 
ness range (Fb cs t — F]im) = 0.08. If any of the 
LMC/SMC proper motion components changes by 
lmasyr -1 , the corresponding model is extremely 
likely to cross the border between the successful 
and insufficient models. If some of the remain- 
ing parameters are altered by further observations, 
our findings concerning the motion of the Clouds 
shall not be affected strongly. 

4.2. Scale Dependence 

What if the volume of the parameter space is 
reduced substantially by changing the proper mo- 
tion ranges of the Clouds? To find the answer, 
let's take a look at Table [TJ It indicates that the 
most influential parameters remain the same, if 
we zoom into the velocity space. The role of the 
choice for the LMC/SMC proper motion is still 
dominant for the evolution of the System. Impact 
of different proper motion components may be al- 
tered if their searched ranges are reduced. It is 
also a notable fact that the values of both func- 
tions Aj and A^p decreased systematically over 
all studied parameters as we switched to the local 
scale of the parameter space. It means that the 
fitness function F does not have a fractal-like (ir- 
regular) structure, and the probability of missing 
steep high peaks (i.e. isolated quality models) in 
the fitness landscape may be reduced by running 
the GA search again on that sub-region of the 
original parameter space, which is of a particular 
interest. We demonstrated that the effective res- 
olution of the GA search may be improved by re- 
ducing the volume of the parameter space. There- 
fore, the exploration of the parameter space for the 
MW-LMC-SMC interaction was performed in two 
levels. First, every LMC/SMC proper motion es- 
timate available was included. Subsequently, the 
proper motio n spread was reduced to the 1-a er- 
ror ranges bv lKallivayalil et all (l2006al lbh to verify 
or eventually correct the global scale search. 



S 



4.3. Spatial motion and tidal models 

We have analyzed the influence of the involved 
parameters on the tidal interaction in the Mag- 
ellanic Syste m. Reproduction of the HI observa- 
tional data ( Briins et al. 2005) by the restricted 
N-body simulation turned out to be critically de- 
pendent on the current spatial velocities of the 
Clouds. 

Since the effective resolution of the search can 
be improved by reducing the volume of the stud- 
ied parameter space (as justified earlier in this 
section), two sets of the LMC/SMC proper mo- 
tion ranges were assumed. After including every 
proper motio n estimate availab l e up to date, only 
the values bv lKallivavalil et al. I (|2006allbT ) were in- 
volved to allow for a high-resolution search on the 
local scale of the parameter space. With the use 
of GA, « 2 • 10 6 parameter combinations, i.e. in- 
dividual N-body simulations, were tested in total, 
and 120 sets providing the highest fitness were col- 
lected for each of the studied volumes of the pa- 
rameter space. 

The resulting models are not considered accept- 
able unless they produce both leading and trailing 
HI Magellanic structures (the Leading Arm and 
the Magellanic Stream). Thus, to allow for quan- 
titative statements, a threshold level of F had to 
be established. In order to do so, one needs to 
understand how the modeled H I morphology and 
kinematics is reflected in the value of the fitness 
function. 

The modeled distribution and kinematics of H I 
in the Leading Arm region and around the main 
LMC and SMC bodies remains similarly unsatis- 
factory over the entire parameter space, especially 
failing to reproduce the observed morphology of 
the Leading Arm (see Fig. [2|). In terms of the fit- 
ness function, the value of F never exceeds 0.3 if 
calculated only for the Leading Arm. 

It is the Magellanic Stream that turned to be 
very sensitive to the choice for the model pa- 
rameters, and critically influencing the resulting 
fitness. Figure [2] illustrates the above discussed 
facts. While the model of F = 0.514 (global 
scale) is able to fit the basic features of the Mag- 
ellanic Stream both in the projected H I distribu- 
tion and the LSR radial velocity profile, the best 
model for the increased LMC/SMC spatial veloci- 
ties (F — 0.336, local scale) places the Magellanic 



Stream to the position-position-velocity space in- 
correctly, as the m — 400 km s -1 tip of the Stream 
is shifted towards the Clouds by ps 20° compared 
to the observations (1 pixel in Fig.[3]equals roughly 
to 2° in the plane of sky), and the modeled mor- 
phology differs seriously from the observed struc- 
ture of the Magellanic Stream. Generally speak- 
ing, the described behavior of the modeled trailing 
stream is responsible for the resulting fitness of a 
given model, and was used to define the desired 
threshold level of the F value. 

Unlike the local scale analysis, the global search 
of the entire parameter space always resulted in a 
model placing the high negative radial velocity tip 
of the Stream to the correct projected position, 
and the modeled trailing stream filled roughly the 
same area of the HI data-cube as the observed 
Magellanic Stream. Following the mentioned fact, 
we selected the worst of the 120 fits from the global 
scale search to represent the threshold level of the 
fitness and so F > 0.434 defines satisfactory mod- 
els of the Magellanic System. 

Figure [3] indicates, that the sub-region of the 
(u N 



w 

l lmc 



)-plane introduced bv lKallivavalil et"al 



(|2006af ) does not allow for any satisfactory models 
assuming the pure tidal interaction. Such a con- 
clusion is also confirmed by the local scale search 
restricted to the HST proper motion data only 
(lower plot of Figure [3]). Thus, Fig. [3] denotes 
that it is not possible to simulate the evolution of 
the Magellanic System by pure tidal models if the 
spa tial velocity of the LMC is as high as predicted 
by iKallivavalil et all ( 2006a ). However, such a 
speculation can only be confirmed if the entire 
parameter space of the interaction is explored, as 
it is based on the behavior of the fitness function 
F in the neighborhood of the selected point. 

Figure |4] summarizes our results. The value of 
every local peak of the fitness function F (i.e. the 
fitness of every model identified by GA) for the 
interaction of the MW-LMC-SMC system is plot- 
ted as function of the fif mc and /x^ c proper mo- 
tion components for the LMC. While the northern 
component /j^ nc allows for satisfactory modeling 
of the Magellanic System over the entire global 
range, the situation is quite different for the LMC 
proper motion in the western direction (/^uL.)- 
The upper right plot of Figure [4] clearly indicates 
that the tidal scheme fails unless the /x^ c values by 
Kallivavalil et all (|2006ah are reduced by « 40 % 
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to reach /ij^ c > — 1.3 mas vr" 1 . iBesla et al.l (|2007l ) 
have shown that the proper motion component 
fj^ ic controls the Galactocentric spatial velocity 
of the LMC. Hence, one can easily see that the 
pure tidal models of the Magellanic System con- 
stantly fail for the spatial velocities putting the 
LMC on highly eccentri c or maybe unbo und or- 
bits about the Galaxy (jBesla et al.l 120071 ) . This 
conclusion was verified and confirmed by the lo- 
cal scale parame ter study focusing on the proper 
motion ranges bv lKallivavalil et al.l(|2006allbT ). De- 
spite our efforts, no even decent models were found 
over this portion of the parameter space (see the 
lower row of Figure |4|) . Please note also that the 
concentration of the GA fits towards lower values 
of fJ^ c still remains even for the reduced velocity 
ranges. 

Within the original volume of the parameter 
space, assuming the tidal interaction, no satisfac- 
tory r eproduction of the H I data by iBriins et al 



(|2005l ) was achieved for the HST proper motions. 
In agreement with the previous studies, the model 
succeeded only if the Clouds were moving at sub- 
stantially lower Galactocentric velocities. The fol- 
lowing high-resolution analysis of the local scale 
of the parameter space (the proper motion ranges 
for the Clouds were restricted to the HST veloc- 
ity data) did not change the previous result and 
no specific quality parameter combinations were 
revealed. Regarding the above summarized facts 
and results, we conclude that the dynamical evo- 
lution of the Magellanic System with the currently 
accepted total mass of the MW cannot be ex- 
plained in the framework of pure tidal models 
for the proper motion data by iKallivayalil et al" 
(|2006al lbh within their 1 —a errors. 

5. Conclusions 



The new results introduced bv lKallivavalil et"al 
(|2006al |bT) have serious consequences for our under- 
standing of the dynamical evolution of the Magel- 
lanic System, and in a wider context they influence 
our view of the Local Group and its formation. 
Such facts, together with the key result of this pa- 
per indicating a conflict between the tidal models 
and observations of the Magellanic System, lead 
us necessarily to the questions about the reliabil- 
ity of the original HST data and correctn ess of 
their treatment by iKallivavalil et al. ([2006, 



Unfortunately, the first issue cannot be ad- 
dressed until the HST measurements are cross- 
checked by an instrument of a competitive resolu- 
tion and other significant physical characteristics. 
In a close future, the GAIA mission should allow 
for high-precision astrometry, and hopefully the 
Magellanic Clouds and their proper motions will 
become objects of its interest as soon as possible. 

Recently, an interesting constraint on the cur- 
rent velocity of the LMC w as introduced by 
McClure-Griffiths et all (j2008l ). They were able 
to estimate the Galactocentric distance to the 
cross-section of the Magellanic Leading Arm with 
the gaseous disk of the MW. Although the ob- 
served part of the Leading Arm does not necessar- 
ily trace the future orbit of the LMC (the Leading 
Arm is believed to lead the Magellanic System), 
the measured position puts the lower limit on the 
Galactocentric distance to the point of the next 
passage o f the LMC through the Galact ic plane. 



Following iMcClure-Griffiths et al.l (|2008l) . the es- 



timated distance is not at odds with the value 
ob tained by adopting the c urrent LMC velocity 
bv IKallivayalil et al. ( 2006a ). In addition to the 
above discussed issue, a verification of the dis- 
tances to the main bodies of the Clouds may be 
worth a consideration. Space velocities also de- 
pend on the distances to the LMC/SMC, which 
should be checked carefully in future to clarify 
how far they may influence our conclusions. 

Concern ing the HST proper moti on data pro- 
cessing by IKallivayalil et al.l ( 2006al lbh. both pa- 
pers indicate that the data were analyzed and in- 
terpreted very carefully and the adopted method 
is well .justifie d. A very strong suppor t to the 

came 



conclusions by IKallivayalil et al. J 2006afbT) 



(2008) 



cur- 



recently from the study by iPiatek et al 
who also derived the LMC and the SMC 
rent proper motion s from the original HST data. 
Piatek et al. ( 20081 ) introduced a different method 



to process the LMC/SMC stellar proper mo- 
tions, b ut their results a gree w ith the find- 
ings by IKallivayalil et al.l ( 2006al lbl) very well. 
They were able to reduce the LMC/SMC proper 
motion errors by the fa ctor of 3 compared to 



Kallivavalil et al.l (|2006al lbl) due to their modi- 
fied treatment of the original data. However, 
while their 1-tr proper motion errors of the LMC 
are completely embed ded within the cor r espond - 



ing error estimates by IKallivavalil et al. ( 2006al ). 
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the mean value of SMC proper motion compo- 
nent /z^L is offset bv_+0.41 masyr -1 compared to 
Kallivavalil et all (|2006bl ) (see Fig.EJ). 

It means that the region of the parameter 
space d elimited by the \-a proper motion er- 
rors by IPiatek et al. ( 2008 ) entered our global 
scale analysis, but it was not explored by the 
local scale search. Nevertheless, there is no in- 
dication that such a high-resolution exploration 
might alter significantly the results obtained for 
the global scale searc h. Since the only differen ce 
between the re sults by Kallivavalil et al.l (|2006al lbl) 
and those by IPiatek et al.1 (|2008l ) exists for the 
western proper motion component of the SMC, 
the Eq. PD|) yields a typical growth in the mod- 
els' fi tness of AF(m^L ) ~ +0. 03 as one switches 
from 



Kallivavalil et al 



(2008), because the /i, 
20081) is bv » 0.4 m a s yr 



w 

sine 



( 2006allbT ) to 

estimate by 



Piatek et al 



Piatek et al 



lower than the result 



bv iKallivavalil et alj (l2006bl). The overall qual 



ity of the models for the Piatek et al.l (|2008l ) data 
is then well below the fitness threshold level es- 
tablished in Sec. 14.31 Under such conditions it is 
extremely unlikely to expect any significant revi- 
sion of the global scale parameter search by the 
high resolution ana lysis of the proper motions by 



Piatek et all (|2008[) 



If we admit that the latest observational esti- 
mates of the current proper motion for both Mag- 
ellanic Clouds are correct, what would be the im- 
pact on our understanding of the Magellanic Sys- 
tem and its evolution? We showed that the pure 
tidal stripping is insufficient to redistribute the H I 
gas in the System to confor m the available ob- 
servations ( Brims et al.ll2005l). Nevertheless, the 
classical model scheme by Fuiimoto fc Sofud ( 19761 ) 
offers an extremely simplified view of the physical 
processes influencing galactic interacti ons. Our re- 
sults, together with the conclusions by lBesla et al 



(|20071 ). do not allow for more than little doubts, 
that despite the high spatial velocities of the 
Clouds, the models of the interaction may succeed 
if sufficiently efficient physical mechanisms are in- 
troduced. 

In the first order, the ram pressure stripping 
(jMastropietro et al.l 120051 ) scenario is well worth 



further efforts, since its efficiency strongly depends 
on several only weakly constrained parameters, 
namely the density of the extended gaseous halo 
of the Galaxy. Moreover, the ram pressure force 



scales as v 2 with the relative velocity of the in- 
teracting objects, an d so the high spat i al velq ci- 
ties of the Clouds by IKallivavalil et all (|2006al lbl) 
may compensate the effect of the corresponding 
reduced timcscalc for the interaction with the 
gaseous halo o f the G alaxy. Recently, the paper by 
Nidever et alJ ( 20081 ) offered an exciting alterna- 
tive to the tidal/ram pressure models, considering 
the possible intense ejection of mass from several 
star-forming regions in the LMC super-shells. 
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Table 1 

Parameter dependence of the fitness function 





Global scale 






Local scale 




3 


Aj ■ 10 2 


Aj,f ■ 10 2 


3 


A, • 10 2 


A jiF ■ 10 2 


N 
Ms™ 


15.96 


31.03 




3.09 


5.53 




8.76 


17.20 




3.01 


5.48 




8.55 


18.92 


Ac 


2.64 


4.85 


N 


6.98 


15.24 


fjL 


2.63 


4.76 


Otoe 


4.22 


8.34 


dime 


1.77 


3.31 


(m - M) amc 


3.95 


8.02 


r->disk 
K lmc 


1.69 


2.94 


(m - M)i mc 


3.92 


7.91 


Mi mc 


1.69 


3.11 


M lmc 


3.86 


7.78 


eimc 


1.67 


3.09 


M smc 


3.48 


6.93 


(m - M) ln 


1.65 


3.06 


A 


3.46 


6.95 


A 


1.63 


2.98 




3.31 


6.52 


(m - M) ar 


1.59 


2.95 


time 


2.60 


5.36 
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Fig. 1. — The 2D projections of the Magellanic parameter space to the (fi , /x )-plane for both the LMC 
(left plot) and the SMC. The gray fillings mark the proper motion ranges explor ed by the GAs. The label s 



indicate t he proper motion s as exp ected by various studie s . K0 6a stands for IKallivavalil et al. ( 2006a[ ). 
K06b for IKallivavalil et all (l2006bh. PI08 for IPiatek et all <l2008l). J94 for Ijones et alj (1 1994). PPM for 



Kroupa et alj (Il994). HIP for |Kro upa fe Bastianl dl997l) P02 for IPedreros i et al.l (|2002[). AKF for the value 
combining Freire et al. ( 2003 ) with I Anderson &; Kind (|2004al ibh . The ellipses show the 68.3% confidence 



regions. 
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Fig. 2. — Visualization of the entire 3D HI data-cube of the Magellanic System. The column density 
isosurface £ = 10~ 5 E max is shown in every plot, together with the data-cubes projected to the 2D maps of 
the integrated column density both in the position-position and position-radial velocity spaces. The best 

—0.34 mas yr 



model ever found in the global scale search (F = 0.514, /iuL. 
hand plot) is compared to the best model identified for the iKallivavalil et al 




—0.70 mas yr -1 , left 
data, i.e. the proper 

motion region delimited by Eq. ([5]) to ([5]), (F = 0.336, nf mc = — l.QS masyr -1 , J = +0 .40masyr~ 1 , right 
hand plot), and to the low-resolution compilation of the HI data bv lBriins et al.l (2005). The dashed line 
in the position-velocity projections depicts the mean radial velocity gradient along the observed Magellanic 
Stream. 
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Fig. 3. — The (fJ.f mc , /x^J-planes for the 2D projections of the function F. For the moment, the remaining 
parameters are fixed to the values corresponding to the best models for the original (upper plot) or the 
reduced proper motion ranges, respectively. The best models are marked by the white c rosses . The white 



ellipse i n the upper plot indicates the 68.3 % confidence region of the LMC proper motion bv lKallivavalil et al 
(|2006al ) and roughly corresponds to the entire {fif mc , ^^ lc )-plane depicted in the lower plot. 
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Fig. 4. — Distribution of all GA fits of the Magellanic System over the analyzed ranges for the LMC proper 
motion components i~if mc and /i^, c . The upper row presents the low-resolution search of the original volume 
of the parameter space. The gray-filled areas indicate the reduced proper motion intervals. They were 
studied subsequently and the resulting span of the 120 identified GA fits is depicted in the lower row. 
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